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Abstract. The metal-rich bulge globular cluster NGC 6388 shows a distinct blue horizontal-branch tail in its colour-magnitude 
diagram (Rich et al. 1997) and is thus a strong case of the well-known 2'"' Parameter Problem. In addition, its horizontal 
branch (HB) shows an upward tilt toward bluer colours, which cannot be explained by canonical evolutionary models. Several 
noncanonical scenarios have been proposed to explain these puzzling observations. In order to test the predictions of these 
scenarios, we have obtained medium resolution spectra to determine the atmospheric parameters of a sample of the blue HB 
stars in NGC 6388. Using the medium resolution spectra, we determine effective temperatures, surface gravities and helium 
abundances by fitting the observed Balmer and helium lines with appropriate theoretical stellar spectra. As we know the distance 
to the cluster, we can verify our results by determining masses for the stars. During the data reduction we took special care 
in subtracting the background, which is dominated by the overlapping spectra of cool stars. The cool blue tail stars in our 
sample with effective temperatures T^ff a 10,000 K have lower than canonical surface gravities, suggesting that these stars are, 
on average, a;0?4 brighter than canonical HB stars in agreement with the observed upward slope of the HB in NGC 6388. 
Moreover, the mean mass of these stars agrees well with theoretical predictions. In contrast, the hot blue tail stars in our 
sample with Teff > 12,000 K show significantly lower surface gravities than predicted by any scenario, which can reproduce the 
photometric observations. Their masses are also too low by about a factor of 2 compared to theoretical predictions. The physical 
parameters of the blue HB stars near 10,000 K support the helium pollution scenario. The low gravities and masses of the hot 
blue tail stars, however, are probably caused by problems with the data reduction, most likely due to remaining background 
light in the spectra, which would affect the fainter hot blue tail stars much more strongly than the brighter cool blue tail stars. 
Our study of the hot blue tail stars in NGC 6388 illustrates the obstacles which are encountered when attempting to determine 
the atmospheric parameters of hot HB stars in very crowded fields using ground-based observations. We discuss these obstacles 
and offer possible solutions for future projects. 

Key words. Stars: horizontal branch - Stars: evolution - Techniques: spectroscopic - Galaxy: bulge - globular clusters: indi- 
vidual: NGC 6388 



1. Introduction 

Ever since its discovery over 30 years ago (Sandage & Wildey 
119671 van den Bergh 1967 1, the 2"*^ parameter effect has stood 
as one of the major unsolved challenges in the study of the 
Galactic globular clusters. While it was recognized quite early 
that the horizontal branch (HB) becomes redder on average 
with increasing metallicity, many pairs of globular clusters are 
known with identical metallicities but markedly different HB 

* Based on observations with the ESO Very Large Telescope at 
Paranal Observatory, Chile (proposal ID 69.D-0231(A)) 
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morphologies, e.g., M 3 versus M 13. Thus some parameter(s) 
besides metallicity (the P' parameter) must affect the evolution 
of the HB stars in these globular clusters. Possible 2"'* parame- 
ter candidates include the globular cluster age, mass loss along 
the red-giant branch (RGB), helium abundance Y, a-element 
abundance, cluster dynamics, stellar rotation, deep mixing, etc. 

Hubble Space Telescope observations by Rich et al. (I1997> 
have found an unexpected population of hot HB stars in the 
metal-rich globular clusters NGC 6388 and NGC 6441 ([Fe/H] 
~ -0.5), making these clusters the most metal-rich clusters to 
show the 2"'^ parameter effect. Ordinarily metal-rich globular 
clusters have only a red HB clump. However, NGC 6388 and 
NGC 6441 possess extended blue HB tails containing »15% 
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of the total HB population. Quite remarkably, the HBs in both 
clusters slope upward with decreasing B-V with the stars at the 
top of the blue tail being nearly (J?5 brighter in V than the well- 
populated red HB clump. Moreover, the RR Lyrae variables in 
these clusters have unusually long periods for the cluster metal- 
licity, leading Pritzl et al. (2000) to suggest that NGC 6388 and 
NGC 6441 may represent a new Oosterhoff group. For all of 
these reasons the HBs of NGC 6388 and NGC 6441 are truly 
exceptional. 

In the next section we will discuss the implications of 
NGC 6388 and NGC 6441 for the 2"'^ parameter effect and 
will review a number of scenarios for explaining the HB mor- 
phology of these clusters. Sect. |3l describes the observations 
and reduction of the medium resolution spectra that we have 
obtained to test these scenarios, while Sect. 0] gives the atmo- 
spheric parameters derived from these spectra. These results 
are compared with the theoretical HB tracks in Sect. ID Sect.|6l 
summarizes our conclusions. 

2. Horizontal-Branch Morphology: Problems and 
Scenarios 

The presence of hot HB stars in globular clusters as metal- 
rich as NGC 6388 and NGC 6441 may provide an important 
diagnostic for understanding the 2"'' parameter effect for the 
following reason. In intermediate-metallicity globular clusters 
such as M 3 the HB spans a wide range in color that extends 
both blueward and redward of the instability strip. The loca- 
tion of a star along the HB is then quite sensitive to changes 
in the stellar parameters. In fact, this is why the HB is "hori- 
zontal". In metal-rich globular clusters, however, the situation 
is different. Due to their high envelope opacity, metal-rich HB 
stars are normally confined to a red clump. To move such stars 
blueward requires a larger change in the stellar structure than 
for intermediate-metallicity stars. Thus any 2""* parameter can- 
didate capable of producing hot HB stars in a metal-rich glob- 
ular cluster might also have other observational consequences. 
Indeed, the upward sloping HBs in NGC 6388 and NGC 6441 
suggest that the 2"^ parameter in these globular clusters is af- 
fecting both the temperature and the luminosity of the HB stars. 

Can canonical models explain the upward sloping HBs in 
NGC 6388 and NGC 6441 ? In principle, one could produce hot 
HB stars in these globular clusters by increasing the cluster age 
or by enhancing the amount of mass loss along the RGB. Rich 
et al. J1997> considered both of these possibilities but found 
neither of them to be satisfactory because the required increase 
in the cluster age is quite large and because the frequency of 
stellar interactions within the cores of these clusters seems too 
low to produce the additional RGB mass loss. This conclu- 
sion was further supported by the theoretical HB simulations 
of Sweigart & Catelan ( 1998 hereafter SC98). They found that 
the HB morphology predicted by canonical HB models is flat 
in the (My, B -V) plane. Increasing the cluster age or the RGB 
mass loss simply moves the models blueward in B -V without 
increasing their luminosity. Thus canonical HB models cannot 
account for the HB morphology of NGC 6388 and NGC 6441. 
In particular, two of the most prominent 2"^^ parameter candi- 
dates - age and RGB mass loss - do not work. 



More recently, Raimondo et al. (12002'* have argued that 
canonical HB models can produce upward sloping HBs if the 
metallicity is sufficiently high or if the mixing-length ratio a 
describing the superadiabatic convection in the outer layers of 
the red HB stars is sufficiently small. Their solar metallicity 
models with a = 1.0 predict a luminosity difference AV be- 
tween the top of the blue tail and the red HB clump that exceeds 
the O'J'5 observed in NGC 6388 and NGC 6441. However, such 
a small value for the mixing-length ratio is untenable because 
it would imply a B - V colour for the RGB that is much red- 
der than observed. The predicted AV for the more appropriate 
mixing-length ratio a = 1 .6 is only x 0^2 for a solar metallicity 
and decreases further with decreasing [Fe/H]. At the metallic- 
ity of NGC 6388 and NGC 6441 the Raimondo et al. models 
predict little difference in the luminosity of HB stars at the top 
of the blue tail compared to those in the red clump, further con- 
firming that some noncanonical effect must be influencing the 
HB morphology in NGC 6388 and NGC 6441. 

The failure of canonical HB models to produce upward 
sloping HBs has prompted the study of other noncanonical so- 
lutions. Theoretical models show that the HB luminosity at a 
fixed metallicity depends on two parameters: the helium abun- 
dance Y and the core mass M^. This fact lead SC98 to suggest 
3 noncanonical scenarios involving increases in either Y or Mc 
which might potentially produce upward sloping HBs. 

The first ("high-F') scenario assumes that the stars in NGC 
6388 and NGC 6441 formed with a high primordial helium 
abundance due to a peculiar chemical enrichment history in 
these clusters. From theoretical models we know that HB tracks 
at high helium abundances have very long blue loops which 
deviate considerably from the zero-age HB (ZAHB). The HB 
simulations of SC98 show that such high-F tracks can in- 
deed produce upward sloping HBs as seen in NGC 6388 and 
NGC 6441 provided Y is very large (;0.4). However, this 
scenario predicts too large a value for the number ratio R of 
HB stars to RGB stars brighter than the HB (Layden et al. 
.1999 ) as well as too bright a luminosity for the RGB bump 
(Raimondo et al. l2002> . Thus a high primordial helium abun- 
dance in NGC 6388 and NGC 6441 can be ruled out. 

The second ("rotation") scenario is based on the fact that in- 
ternal rotation within an RGB star can delay the helium flash, 
thereby leading to a larger core mass and greater mass loss 
near the tip of the RGB. This increase in the core mass Mc 
together with the corresponding decrease in the total mass M 
will shift a star's HB location towards higher effective temper- 
atures and luminosities. HB simulations show that this scenario 
can also produce upward sloping HBs similar to those observed 
in NGC 6388 and NGC 6441. However, it is difficult to un- 
derstand how the blue HB stars in NGC 6388 and NGC 6441 
could have the high rotation rates required by this scenario. 
Preliminary results from high resolution spectra of three cool 
blue HB stars (Teff < 10,000 K) in NGC 6388 do not show ev- 
idence for rotation velocities of v sin / > 10 km sec ' (Moehler 
& Sweigart "20061. Thus the rotation scenario seems unlikely. 

The third ("helium-mixing") scenario is motivated by the 
large star-to-star abundance variations which are found among 
the red-giant stars within individual globular clusters and which 
are sometimes attributed to the mixing of nuclearly processed 
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material from the vicinity of the hydrogen shell out to the stellar 
surface (Kraft 1994 1. The observed enhancements in Al are par- 
ticularly important because they would require the mixing to 
penetrate deeply into the hydrogen shell (Cavallo et al. .l998j . 
Such mixing would dredge up fresh helium together with Al, 
thereby increasing the envelope helium abundance and lead- 
ing to a brighter RGB tip luminosity and hence greater mass 
loss. Thus a helium-mixed star would arrive on the HB with 
both a higher envelope helium abundance and a lower mass 
and would therefore be both bluer and brighter than its canon- 
ical counterpart - just what is needed to produce an upward 
sloping HB. Indeed, the HB simulations of SC98 confirm that 
helium mixing can produce HB morphologies similar to those 
in NGC 6388 and NGC 6441. However, the existence of he- 
lium mixing on the RGB can be questioned on several grounds. 
The O-Na and Mg-Al anticorrelations observed in turnoff stars 
of NGC 6752 by Gratton et al. ( 1200 It indicate that the Al en- 
hancements are more likely due to primordial pollution from an 
earlier generation of stars than to deep mixing on the RGB. It is 
also questionable whether the mixing currents could overcome 
the large gradient in the mean molecular weight within the hy- 
drogen shell of a RGB star. Thus helium mixing also seems 
unlikely. 

These difficulties have prompted a number of additional 
suggestions for explaining the HB morphologies of NGC 6388 
and NGC 6441. One of the earliest suggestions, offered by 
Piotto et al. ( 1997 1, was a spread in metallicity. In this case the 
blue HB stars would be metal-poor compared to the stars in the 
red HB clump. Because the HB becomes brighter in Y with de- 
creasing metallicity, the blue HB stars would also be brighter. 
Thus a metallicity spread might also produce an upward slop- 
ing HB. This possibility was studied by Sweigart (I2002t. who 
showed that upward sloping HBs similar to those in NGC 6388 
and NGC 6441 would require the stars at the top of the blue HB 
tail to be approximately 2 dex more metal-poor than the stars 
in the red HB clump. However, Raimondo et al. (2002 1 have 
noted that the progenitors of the blue HB stars should appear 
as a population of metal-poor giants lying well to the blue of the 
metal-rich RGB. Since such metal-poor giants are not seen in 
the colour-magnitude diagrams of NGC 6388 and NGC 6441, 
Raimondo et al. (I2002> conclude that any metallicity spread 
must be small. Again preliminary results of the high resolution 
spectra mentioned above do not show any evidence for a signif- 
icant metal deficiency in these blue HB stars compared to the 
overall metalHcity of NGC 6388 (Moehler & Sweigart '2006'l. 

Piotto et al. (1997) also suggested that NGC 6388 and 
NGC 6441 might contain two stellar populations with differ- 
ent ages. This possibility has been further explored by Ree et al. 
(2002 1. Their population models for NGC 6388 and NGC 6441 
are able to produce blue HB stars provided these stars are older 
by 1.2 Gyr and metal-poor by 0.15 dex compared to the stars 
in the red HB clump. Such a small difference in metallicity 
between the blue and red HB stars avoids the problem with 
the missing metal-poor giants discussed above. However, while 
such models might produce bimodal HBs, they do not produce 
an upward sloping HB. As shown by SC98, differences in age 
merely move an HB star horizontally along the HB, and, as 
shown by Sweigart (i2002i) . a metallicity difference of 0.15 dex 



is too small to produce a significant HB slope. Thus the Ree 
et al. (2002 1 models do not account for a key property of the 
HBs in NGC 6388 and NGC 6441. Ree et al. ( |200^ also sug- 
gest that the long RR Lyrae periods might be explained if these 
stars are highly evolved from the blue HB. While such stars 
would be brighter and hence have longer periods, they would 
also evolve rapidly across the instability strip on their way to 
the asymptotic-giant branch (AGB). Explaining the observed 
number of RR Lyrae stars under such a scenario would there- 
fore be very difficult (Pritzl et al. l2002t . 
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Fig. 1. The colour-magnitude diagram of the hot horizontal- 
branch stars in NGC 6388 (Piotto et al. ll997> with our spec- 
troscopic targets marked by open circles. 

A more promising possibility is based on the suggestion 
by D'Antona & Caloi (2004) that the stars in globular clusters 
with blue HB tails are born in two events: a first generation 
of helium-normal stars and a second generation of helium-rich 
stars which form from the ejecta of the first generation. Various 
candidates have been proposed for producing this helium-rich 
ejecta, including AGB stars, massive main-sequence stars, type 
II supernovae, and even stars external to the globular cluster 
(Bekki & Norris l2006l and references therein). Support for such 
helium pollution comes from the recent discoveries of a dou- 
ble main sequence in u) Cen (Anderson ,1997. , Norris 2004^ 
Piotto et al. I2005> and a blueward extension of the main se- 
quence in NGC 2808 (D'Antona et al. l2005t . These globular 
clusters apparently contain a significant population of helium- 
rich (Y w 0.4) stars. Since a helium-rich star has a lower turnoff 
mass at a given age, it will be bluer on the HB than a helium- 
normal star. It will also be brighter due to the increased energy 
output of the hydrogen-burning shell. Thus the spread in the 
helium abundance predicted by this helium pollution scenario 
will lead to a spread in color along the HB, with the red clump 
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Table 1. Coordinates and photometric data for our target stars (from Piotto et al. ll997> 



number 




l52000 


V 


B-V 


setup 


1283 


17''36"'23-899 


-44°44'58'.'45 


18"690±0"032 


+0"230±0"'041 


1 


1775 


17h36" 14-921 


-44°43'08'.'87 


18"8 11 ±0^^040 


+0l^92±0"'047 


1 


2612 


17h36"'22-560 


-44°44'01'.'93 


19'^622±0'!'059 


+0'^79±0"'068 


1 


4113 


17h36"'22-884 


-44°44'24'.'63 


16'^76±0"026 


+0"545±0"037 


1 


6239 


n^se"" 19-813 


-44°43'05'.'88 


19"147±0'1'036 


+0'^244±0"'041 


1 


7104 


n^se"" 17-879 


-44°45'05'.'02 


19l'310±0"026 


+0l^30±0"'037 


1 


7788 


17h36"'21-842 


-44°42'55'.'70 


17'1'410±0"023 


■+-0'^360±0"'026 


1 


396 


17''36™21-160 


-44° 44' 24'.' 42 


17'1'383±0"025 


+0l'303±0"te9 


2 


1233 


17''36™24-565 


-44°43'15'.'21 


16'l'867+0"029 


■+-0l'428±0"'038 


2 


2483 


17''36™21-915 


-44°44'08'.'02 


16"927+0"024 


■+-0'^467±0"'028 


2 


2897 


17''36™23-334 


-44° 43' 56'.' 46 


18™835+0™029 


•t-0™264±0"'038 


2 


5235 


17''36™21-301 


-44°45'14'.'05 


17™113+0™019 


•f0™341±0"'024 


2 


7567 


17''36™ 16-614 


-44°44'59'.'74 


18'!'882+0™145 


+0™343±0"'158 


2 



Table 2. Observation parameters for the spectroscopic data. The FORS2 value of the seeing has been measured on the acquisition 
image taken immediately before the spectroscopic exposures. Between the first and second observation of setup 2 no acquisition 
image was taken. 



setup 


date 


start of 


exposure 


seeing 


airmass 


moon 








exposure 


time 


DIMM 


FORS2 




illumination 


distance 


1 


2002-07-10 


03:36:27.049 UT 


3320 sec 


a'9i 


a' 75 


1.095 


0.1% 


153-7 




2002-08-02 


03:37:04.339 UT 


3238 sec 


1'.'78 


1'.'25 


1.245 


42.5% 


136-6 


2 


2002-08-02 


01:34:19.337 UT 


3238 sec 


1'.'70 


I'.'OO 


1.075 


43.5% 


135-5 






02:31:37.199 UT 


3238 sec 


1'.'75 




1.121 


43.0% 


136-0 



stars corresponding to the helium-normal, first generation stars 
and the blue tail stars being progressively more helium-rich as 
the effective temperature increases. An upward sloping HB is a 
natural consequence of this spread in heUum. The fact that the 
HB slope is more prominent in NGC 6388 and NGC 6441 than 
in other blue tail globular clusters may be simply due to their 
higher metallicity which requires a larger increase in helium in 
order to force a star blueward of the red clump. We emphasize 
that this scenario differs from the high-F scenario, mentioned 
above, in which all of the stars are helium-rich and from the 
helium-mixing scenario in which the spread in helium arises 
from deep mixing on the RGB. 



3. Observations and Data Reduction 

3.1. Target Selection 

Our spectroscopic targets were selected from the catalog of 
Piotto et al. J 19971 Fig[l), and include six stars in the roughly 
horizontal part of the horizontal branch (16"6 < V < 17™6) 
and seven stars at the hot end of the blue tail (18™6 < V < 
19™8). We tried to select the most isolated stars from the orig- 
inal WFPC2 images, which were obtained close to the core of 
the globular cluster The coordinates and photometry for our 
targets are given in Tabled Fig.Elshows finding charts for the 
stars. 



All of the scenarios which can produce an upward sloping 
HB predict that the gravities of the blue HB stars should be 
lower than the gravities of canonical blue HB stars. In 1998 we 
observed some of the brighter objects in both clusters (B <18; 
four stars in NGC 6388, three in NGC 6441) to test this pre- 
diction. Unfortunately the results only added to the confusion 
as most of the stars had even higher gravities than predicted 
by canonical evolution (Moehler et al. .1999) . In 2002 we ob- 
tained additional medium resolution spectra of about a dozen 
blue HB stars along the blue tail of NGC 6388 to determine 
their effective temperatures and surface gravities. In the follow- 
ing sections we will describe our analysis of these spectra and 
will compare the derived effective temperatures and gravities 
with the predictions of the helium pollution scenario. 



3.2. Spectroscopy 

We obtained medium-resolution spectra (R sil200) at the VLT- 
UT4 (Yepun) with FORS2 in service mode (see Table |2l for 
details). We used the multi-object spectroscopy (MOS) mode 
of FORS2 (slit length 20"respectively 22") with the standard 
collimator (O'.'2/pixel), a slit width of 0'.'6 and grism B600. The 
slitlets were positioned to cover at least the wavelength range of 
3700 A to 5200 A. As F0RS2 is equipped with an atmospheric 
dispersion corrector, MOS observations at higher airmass are 
not a problem. Due to a mistake during the Phase 2 preparations 
the mask for the second observation of setup 1 was prepared 
with r.'2 wide slitlets. In combination with the mediocre seeing 
these wide slits made the data useless due to the high crowding 
evident in Fig.|2l Therefore we have only one spectrum for each 
star in setup 1 . 



Moehler & Sweigart: Hot HB Stars in NGC 6388 



5 



Fig. 2. Finding charts for all target stars from a FORS2 image 
obtained in white light with a seeing of 0'.'6. 



FORS2 has two 2kx4k thinned MIT CCD detectors with 
anti-reflection coating and a pixel size of (15/vm)^ with a gain 
of 0.7 e /count and a read-out noise of 2.7 e . The CCDs are 
read out with binning 2x2, resulting in a scale of O'.'25/pixel 
respectively 1.5 A/pixel. Due to the limited field-of-view of the 
WFPC2 photometry, from which we selected our targets, we 
used only the master CCD for our observations. 

For each night screen flat fields with two different illu- 
mination patterns and CdHeHg wavelength calibration spec- 
tra were observed. As part of the standard calibration we were 
also provided with masterbias frames for our data. The master- 
bias showed no evidence for hot pixels and was smoothed with 
a 30x30 box filter to keep any possible large scale variations 
while erasing noise. The flat fields were averaged for each night 
and bias-corrected by subtracting the smoothed masterbias of 
that night. From the flat fields we determined the limits of the 
slitlets in spatial direction. Each slitlet was extracted and from 
there on treated like a long-slit spectrum. The flat fields were 
normalized with 4* to 7* order polynomials. The dispersion 
relation was obtained from the wavelength calibration frames 
by fitting 4"* to 6* order polynomials to the line positions along 
the dispersion axis. We used 15 to 19 unblended lines between 
3600 A and 5900 A and achieved a typical rm.s. error of 0.05 A 
to 0.08 A per CCD row. 

Due to exposure times of more than 50 minutes the sci- 
entific observations contained a large number of cosmic ray 
hits. Those were corrected with the algorithm described in Pych 
J2004> with the default parameters. To ensure that this proce- 
dure did not introduce any artefacts, we also reduced the un- 
corrected frames for comparison. We did not find any prob- 
lems with the cosmic ray rejection routine. The slitlets with 
the stellar spectra were extracted in the same way as the flat 
field and wavelength calibration slitlets. The smoothed master- 
bias was subtracted, and the spectra were divided by the corre- 
sponding normalized flat fields, before they were rebinned 2- 
dimensionally to constant wavelength steps. We then corrected 
the curvature of the spectra along the spatial axis using a rou- 
tine by O. Stahl (priv. comm.), which traces a spectrum in a 
predefined spatial region and determines and corrects the cur- 
vature derived that way for the complete slitlet. If possible (i.e. 
if the target was sufficiently isolated and/or sufficiently bright) 
we used the target spectrum for the curve tracing, otherwise the 
brightest spectrum in the slitlet. 

3.3. Sky Subtraction 

The MOS observations were obtained in very crowded fields, 
and due to the fact that several stars were observed simulta- 
neously we could not orient the individual slitlets in a way to 
avoid nearby stars. Therefore most slitlets contain spectra of 
several stars, in many cases overlapping so strongly that it is 
impossible to directly extract the spectrum of our intended tar- 
get (see Fig. |3}. In order to account for this overlap, we pro- 
ceeded as follows: 



1 . we corrected the curvature of the FORS2 spectra 

2. we averaged the wavelength-calibrated two-dimensional 
spectra along their dispersion axis between 3500 A and 
5200 A (roughly the range which is later used for fitting the 
line profiles), thereby producing a one-dimensional spatial 
profile along the slitlet (cf. Fig.|3l upper histogram) 

3. the one dimensional spatial distribution of light was fitted 
with a combination of Moff'at functions, i.e. 

n 

with 



For each profile the parameters aj (amplitude) and bj (posi- 
tion) were fitted individually, whereas the parameters c and 
d, which determine the profile width and shape (and should 
depend only on the seeing and instrumental broadening), 
had to be the same for all objects within one slitlet. We 
used up to 13 individual profiles to fit the full spatial light 
distribution along one slitlet (see Fig.|3j light solid line). 

4. After achieving a good fit to the observed spatial profiles we 
kept the parameters bj, c, and d fixed and used these pro- 
files to fit the spatial profile now at every wavelength step. 
The amplitudes cij and the spatial constant bck (for the true 
sky background) were allowed to vary with wavelength, in 
order to describe the various spectra. 

5. The sum of all profiles except the target profile was then 
used as background (cf. Fig.|3 dashed line) and subtracted 
from the wavelength calibrated two-dimensional image. 
The resulting image - containing only the target spectrum 
- was again averaged over the same wavelength range as 
above to verify the quality of the sky correction (see Fig.|3j 
lower histogram). 

The sky-subtracted spectra were extracted using Home's 
J1986t algorithm as implemented in MIDAS. As is well- 
known, optimum extraction procedures fail for data with good 
S/N. We therefore also extracted the target spectra by simple 
averaging and compared the results from the two procedures 
in order to avoid artefacts introduced by optimum extraction. 
Finally the spectra were corrected for atmospheric extinction 
using the extinction coefficients for La Silla (Tiig l 19771 as im- 
plemented in MIDAS, because they provide the closest approx- 
imation to Paranal conditions, for which no spectroscopic ex- 
tinction coefficients are available. 

For a relative flux calibration we used response curves de- 
rived from spectra of LTT 6248 with the data of Hamuy et al. 
J1992t . The response curves were fit by splines and averaged 
for the two nights. The individual target spectra were corrected 
for Doppler shifts determined from Balmer and helium absorp- 
tion lines and co-added if both available spectra were of com- 
parable quality. Finally all spectra were corrected for an inter- 
stellar extinction of 0™37. In order to allow easy comparison 
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of the spectra, we normalized them using the model spectra 
provided by the fit procedure (see Sect.l^. These normalized 
spectra are shown in Fig.0] 



1000 



500 




4000 



2000 



5 10 10 15 

position along slit ["] 

Fig. 3. The spatial light distribution along two slitlets (left: star 
7104, setup 1; right: star 2483, setup 2). The upper histogram 
is the observed light distribution, the grey continuous line gives 
the fit of all sources, the grey dashed line marks the subtracted 
sky background {- fit of all sources except target). The lower 
histogram is the spatial profile of the sky subtracted image. 

Star 7567 shows a prominent G-band and a clear Mg lb 
absorption, which in combination with its very red colour and 
high photometric error, suggest that this star has a cool com- 
panion (either a physical one or an unresolvable blend). It is 
therefore excluded from the further analysis and discussion. 



4. Analysis 

4.1. Atmospheric Parameters 

For the analysis of the FORS2 spectra we used two sets of 
ATLAS9 model atmospheres (Kurucz I1993> : For stars with 
V < 17™5, which are presumably cooler than 11,000 K and 
therefore not affected by diffusion, we used a metallicity of 
[M/H] = -0.5 and solar helium abundance. For the fainter stars 
withy > 18 "5, which are presumably hotter than 12,000 K and 
therefore affected by diffusion, we used a super-solar metallic- 
ity of [M/H] - +0.5 and helium abundances ranging from solar 
to 1/100 solar. These models should account in a rough way 
for the effects of radiative levitation of heavy elements and 
gravitational settUng of helium (see Moehler et al. 120001 for 
details). From these model atmospheres we calculated spec- 
tra with Lemke's version' of the LINFOR program (devel- 
oped originally by Holweger, Steffen, and Steenbock at Kiel 
University). To establish the best fit, we used the routines de- 
veloped by Bergeron et al. ( 1992.) and Saffer et al. ( 1994) , as 
modified by Napiwotzki et al. (I1999> . which employ a test. 
The cr necessary for the calculation ofx^ is estimated from the 
noise in the continuum regions of the spectra. The fit program 
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Fig. 4. Normalized spectra for all target stars. Spectral regions 
affected by cosmics are left out. On the right we show the re- 
gion of the Mg lb triplet, which is quite prominent in cool stars 
and may indicate contamination of our spectra by such stars. 



normalizes model spectra and observed spectra using the same 
points for the continuum definition. We used the Balmer lines 
Hyj to Hi2 (excluding to avoid the Ca n H line) for the fit. For 
the hot stars we included also the He i lines AA 4026 A, 4388 A, 
4471 A, 4921 A to determine helium abundances, whereas the 
helium abundance was kept fixed at the solar value for the cool 
stars. The results are given in Tableland plotted in Fig.|5l 

Recent tests have shown, however, that these fit routines 
underestimate the formal errors by at least a factor of 2 
(Napiwotzki 2005, priv. comm.). In addition, the errors pro- 
vided by the fit routine do not include possible systematic er- 
rors due to, e.g., flat field inaccuracies or imperfect sky subtrac- 
tion. As we cannot estimate these systematic errors, we multi- 
plied the formal errors by 2 to provide at least a lower limit to 
the errors. 



4.2. Reddening 



Using the effective temperatures derived this way, we also es- 
timated reddenings by comparing the predicted colour from 
' For a description see,http://a400.sternwarte.uni-erlangen.de/~ai26/linfli^llnfcMLBlint del colours (for the same metaUicities as the model 
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Table 3. The atmospheric parameters and estimated reddenings 
from the spectroscopic analyses of our target stars. The spec- 
troscopic results for stars below 9000 K (italized in table) are 
considered uncertain. For all stars below 11, 000 K we assumed 
a solar atmospheric helium abundance, as these stars should not 
be affected by diffusion (see text for details). 



star 


Teff 


logg 


log^ 


pspec 

B-V 




[K] 


[cm s"^] 




396 


10200+ 300 


3.52+0.12 




0™360 


1233 


7640± 110 


2.85+0.72 




0'"2S7 


1283 


17200+ 760 


4.28+0.12 


-1.71+0.22 


0™408 


1775 


16100+ 880 


4.04+0.18 


-1.68+0.32 


0";'464 


2483 


7780± 110 


2.97±0.12 




0^^345 


2612 


24800+3300 


4.68+0.38 


-1.66+0.20 


0";"511 


2897 


18700+1400 


4.59+0.20 


-1.48+0.32 


0"i"451 


4113 


7200± 140 


2.72±0.22 






5235 


9960+ 320 


3.61+0.16 




0^385 


6239 


19300+1900 


4.33+0.26 


-1.55+0.26 


0"443 


7104 


19500+2000 


4.43+0.24 


-2.06+0.36 


0™429 


7788 


10200+ 280 


3.51+0.14 




0"360 



spectra used in the fits) to the observed colours. Table|3shows 
an obvious and very suspicious trend, namely, that the redden- 
ing increases towards higher effective temperatures: The cool 
stars show an average reddening of 0™35+0™05, which is con- 
sistent with the literature value for NGC 6388. However, the 
fit results for stars below 9000 K are uncertain, as such stars 
show more and more metal lines, which are not included in the 
model spectra and possibly not evident in the observed spec- 
tra due to the low resolution of the data (although some lines 
can be seen in the spectra of the reddest stars 4113, 2483, and 
1233 in Fig.|3}. Restricting the averaging to the stars around 
10,000 K yields a mean reddening of 0'"39. The hot stars, how- 
ever, show an average reddening of 0"45+0"04, which is con- 
siderably larger than the accepted value. While Raimondo et al. 
J2002> report on possible differential reddening of up to 0™1, 
the correlation between reddening and temperature looks sus- 
picious and might indicate an overestimate of the effective tem- 
peratures in the hot stars. As discussed below such an overesti- 
mate may be due to a partial filling of the line cores by residual 
light from cool stars (see Sect.|^). 

At least a part of the correlation might be due to the differ- 
ent effective temperatures of the stars. As discussed by Grebel 
& Roberts (1995 ) and others, the reddening effect of interstel- 
lar matter depends on the spectral type of a star. On average the 
reddening increases for bluer spectral types. Thus one would 
expect blue HB stars to have a higher reddening than the red 
giants in the same globular cluster However, using their Fig. 9, 
we estimate that the reddening difference between red giants 
with [M/H] - -0.5 and blue tail stars with an atmospheric 
metallicity of [M/H] - +0.5 is at most about 0"04 for a red- 
dening of Eb-v - 0"37 and the difference between cool blue 
HB stars and hot blue tail stars is below 0™01. This effect could 
therefore explain the slightly higher reddening of the cooler 
blue HB stars compared to the reddening of the cluster in gen- 
eral. 



However, from the colour-magnitude diagram in Fig.^ it is 
obvious that our "hot" targets have a higher reddening than the 
"cool" ones, as they lie all on the red side of the blue tail and 
in some cases even overlap with the cooler and brighter targets 
in B-V. 



5. Comparison with Theoretical Tracks 




80000 10000 8000 

■I'c, [K] 

Fig. 5. Effective temperatures and surface gravities for our tar- 
get stars as derived from line profile fits. For comparison we 
show the zero-age HB (ZAHB) for both a canonical helium 
abundance (7 = 0.23, solid line) and two helium-rich composi- 
tions {Y - 0.33, long dashed line; Y - 0.43, short dashed line). 
For comparison we also show the results for the metal-poor 
globular cluster NGC 6752 from Moehler et al. ( l2000l . 

We next compare the derived atmospheric parameters given 
in Table |3] with the theoretical HB tracks for both canonical 
and helium-rich compositions. Our goals are to determine if the 
blue tail stars provide any support for the helium enhancement 
predicted by the helium pollution scenario and, more specifi- 
cally, to see if the stars near the top of the blue tail have the 
lower than canonical gravities inferred from the upward slope 
of the HB in NGC 6388. In Fig.|5]we plot the effective temper- 
atures and surface gravities of our target stars together with the 
ZAHBs for a canonical helium abundance (Y - 0.23) and two 
helium-rich compositions {Y - 0.33, 0.43), as might arise from 
helium pollution. According to the helium pollution scenario, 
the helium abundance should increase along the HB from its 
canonical value in the reddest HB stars to a helium-rich value 
along the blue tail. Thus the blue tail stars should show the 
maximum offset from the canonical surface gravities. 

The stars in Fig.|5]fall into three groups: three stars cooler 
than 9000 K, three stars near 10,000 K and six stars hotter than 
12,000 K. As discussed in the previous section, the coolest stars 
with Teff below 9000 K have less reliable parameters and thus 
will not be considered further. 
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5.1. Cool Blue Tail Stars Near 10,000 K 

The moderately cool stars near 10,000 K lie above the canon- 
ical ZAHB in Fig. |5] At first glance this might be interpreted 
as evidence for a helium enhancement in these stars of per- 
haps as much as Y x 0.33. However, before drawing such a 
conclusion, we need to determine if the masses of these stars, 
as obtained from their atmospheric parameters, are consistent 
with the theoretical HB masses. Previous experience with the 
analysis of hot stars in globular clusters has shown that low sur- 
face gravities are often associated with impossibly small HB 
masses, thus indicating that the surface gravities are not reli- 
able. To carry out this consistency check, we have estimated 
masses for all of our target stars using the approach outlined by 
Moehler et al. (2000l. Fig. |6l compares these masses with the 
theoretical ZAHB masses for both canonical and helium-rich 
compositions. We note that the ZAHB masses in Fig.|6ldepend 
only weakly on the helium abundance. The mean mass of the 
cool blue tail stars is 0.52 in good agreement with the the- 
oretical masses. 

Given this consistency check, we can now examine how 
the predicted surface gravity and luminosity of the cool blue 
tail stars depend on the helium abundance. In order to reduce 
the uncertainty in our results, we will use the average effective 
temperature and surface gravity of these stars instead of the 
individual measurements, i.e., <Teif > = 10,120 K and <log g> 
= 3.547. These average values should better represent the mean 
atmospheric parameters during the HB phase. The dependence 
of the surface gravity at Tefr - 10,120 K on F is given in the 
top panel of Fig.0for stars at the ZAHB, at the midpoint of the 
HB phase and at the point 90% through the HB phase. We see 
that the bulk of the HB phase is spent within a narrow range 
in log g at a given helium abundance. Only near the end of the 
HB phase does the surface gravity decrease significantly below 
its ZAHB value. The corresponding variation in the luminosity 
with Y is shown in the lower panel of Fig.0 We again see that 
most of the HB phase is spent close to the ZAHB luminosity 
as long as F < 0.35. At higher helium abundances an HB star 
will evolve along a blue loop that becomes increasingly more 
extended in effective temperature and that deviates more and 
more from the ZAHB. This leads to the substantial increase in 
the predicted luminosity width of the HB evident in Fig.0 

The results in Fig.0are combined in Fig.^to show the pre- 
dicted variation of the luminosity with surface gravity at Tgff = 
10,120 K. The thin vertical line represents the mean surface 
gravity of our cool blue tail stars. We see that the luminosity 
at this surface gravity does not depend very much on whether 
the cool blue tail stars are near the ZAHB or near the end of 
the HB phase. Assuming that, on average, they are near the 
midpoint of the HB phase, we find a mean luminosity <log 
L> of 1.62 for log g - 3.547. The corresponding luminosity of 
a canonical HB star at the midpoint of its HB phase is log L 
= 1.46 at Tetf - 10,120 K. The difference in these luminosi- 
ties implies that our cool blue tail stars are ^0^4 brighter than 
canonical HB stars. This result is in reasonable agreement with 
the (J^5 increase in luminosity implied by the upward slope of 
the HB in NGC 6388. While our sample of cool blue tail stars 
is admittedly very small and the errors in their surface gravities 



significant, our results do suggest a brighter than canonical lu- 
minosity that is at least consistent with the scenarios that can 
explain the HB slope in NGC 6388. 

One can also ask what helium abundance is required to pro- 
duce the lower surface gravities of the cool blue tail stars. From 
the top panel of Fig. Qwe find that a mean surface gravity of 
log g = 3.547 corresponds to a helium abundance of F = 0.32 
if the cool blue tail stars are in the main part of their HB phase. 
If, however, they are near the end of their HB phase, then the 
required helium abundance would be F = 0.28. In either case 
a helium enrichment, as might be expected from the helium 
pollution scenario, could account for the properties of the cool 
blue tail stars in our sample. Accounting for such an increased 
helium abundance in the theoretical model spectra, which were 
used to determine effective temperatures and surface gravities, 
would lower the derived surface gravities, because an increase 
in helium abundance broadens the Balmer line profiles. Thus a 
given observed profile will yield a lower surface gravity if fitted 
with helium-enriched model spectra. 
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Fig. 6. Effective temperatures and masses for our target stars. 
For comparison we show both the canonical ZAHB (F = 0.23, 
solid line) and two helium-rich ZAHBs (F - 0.33, long dashed 
line; F - 0.43, short dashed Une). 



5.2. Blue Tail Stars Hotter Than 12,000 K 

The hot blue tail stars in Fig. [S] have surface gravities that are 
substantially smaller than the canonical values and indeed even 
smaller than the surface gravities along the helium-rich ZAHB 
for F - 0.43. As in the case of the cool blue tail stars, we can 
use the masses obtained from the atmospheric parameters to 
test the reliability of these surface gravities. The masses of the 
hot blue tail stars plotted in Fig.|6]are, on average, about a fac- 
tor of 2 smaller than the theoretically predicted masses and, 
moreover, are strongly correlated with the offset in log g from 
the ZAHB in Fig.|5l The only exception is star 2897 at Tefr = 
18,700 K, which has a mass in good agreement with the the- 
oretical value and which lies close to the F - 0.43 ZAHB in 
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Fig. 7. Predicted surface gravity and luminosity of HB stars at 
TefF = 10,120 K as a function of the helium abundance Y. The 
solid curves refer to the ZAHB, the dashed curves to the mid- 
point of the HB phase, and the dotted curves to the point 90% 
through the HB phase. 



Fig. Is] The atmospheric parameters of star 2897, if representa- 
tive of the true atmospheric parameters of a typical hot blue tail 
star, would be consistent with an increase in the helium abun- 
dance along the blue tail. However, the surface gravities of the 
other hot blue tail stars seem spuriously too low. 

The fact that the cool blue tail stars, which are visually 
about 2 magnitudes brighter than the hot blue tail stars, have 
reasonable parameters in Figs. 5 and 6 suggests that the prob- 
lems with the parameters of the hot blue tail stars may be due to 
problems with background subtraction. A possible effect would 
be a filling of the Balmer line cores by residual light from cool 
stars, which would simulate a higher effective temperature (= 
less deep cores) for the hotter stars. The wings would be less 
affected by such light, and the surface gravity determined from 
them would probably be closer to the true value. However, with 
our current data we have no way of verifying this possibility. 

For comparison we show also results for stars in the metal- 
poor globular cluster NGC 6752, which were analysed in the 
same way as the stars discussed here, but are almost un- 
crowded. Obviously the stars in NGC 6388 show lower sur- 
face gravities at a given temperature than those in NGC 6752, 
whereas canonical stellar evolutionary theory would argue in 
the opposite direction. This may be taken as support for the in- 
fluence of background subtraction. On the other hand, however, 
one should keep in mind that the helium enrichment required 
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Fig. 8. Variation of the luminosity of HB stars at Tefi = 10,120 
K with surface gravity over a range in the helium abundance 
from Y = 0.23 to F = 0.43. The soUd curve refers to the ZAHB, 
the dashed curve to the midpoint of the HB phase, and the dot- 
ted curve to the point 90% through the HB phase. The values 
of Y, indicated by the solid squares along the ZAHB curve, in- 
crease from Y = 0.23 in the lower right to F = 0.43 in the upper 
left in increments of 0.05. The thin vertical line at log g - 3.547 
denotes the mean gravity of the cool blue tail stars. 



to produce hot HB stars in a metal-rich globular cluster would 
produce significantly higher luminosities and thus lower gravi- 
ties. 

Blue tail stars hotter than 12,000 K are also affected by ra- 
diative levitation which can increase the atmospheric metallic- 
ity to super-solar values (Behr et al. ll999l Moehler et al. 120001 
Moehler l2001l Behr '20031. We attempted to include this effect 
in our analysis by using metal-rich stellar atmospheres when 
determining the atmospheric parameters of the hot blue tail 
stars. However, these stellar atmospheres assume a scaled so- 
lar abundance distribution which may not adequately account 
for the highly non-solar abundance distribution produced by 
radiative levitation. Further study is needed to see if this differ- 
ence between the assumed and actual abundance distributions 
can partly explain the low gravities of the hot blue tail stars. 
This explanation is supported by the low gravities found for 
hot blue tail stars in NGC 6752 (Moehler et al. 20Q0j and M 13 
(Moehler et al. l2003> . It is important to note that the low he- 
lium abundances observed in the hot blue tail stars are caused 
by diffusion effects and do not contradict any helium pollution 
scenario. 

Fig. |5] compares the surface gravities of the hot blue tail 
stars only with the surface gravities of ZAHB models. Since at 
high helium abundances the HB evolutionary tracks can deviate 
considerably from the ZAHB, we also need to examine how the 
post-ZAHB evolution affects the predicted gravities. In Fig.|9] 
we overplot the hot blue tail stars onto an HB simulation cover- 
ing the entire HB phase for Y - 0.43. For simplicity a uniform 
distribution in mass was assumed in this simulation. The hottest 
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Star in our sample (star 2612) is not included in Fig.|5]because 
its surface gravity is much too low to be explained by HB evo- 
lution. Not surprisingly, star 2897 lies within the bulk of the 
HB stars in Fig.|5] However, all of the other hot blue tail stars 
lie in regions that are poorly populated. We conclude that HB 
evolution is unlikely to account for the low surface gravities of 
the hot blue tail stars - a further indication of the uncertainty in 
the atmospheric parameters of these stars. 



3.8 




22000 20000 18000 16000 
Teff 

Fig. 9. Comparison of the atmospheric parameters of the hot 
blue tail stars with an HB simulation for a helium abundance Y 
= 0.43. 



6. Conclusions 

Our results for the cool blue tail stars in NGC 6388 are con- 
sistent with both the upward slope of the HB in NGC 6388 
and the predictions of the helium pollution scenario. Moreover, 
they resolve the puzzling conundrum discussed by Moehler et 
al. ( I1999i who found higher than canonical surface gravities 
for a sample of blue tail stars in NGC 6388 and NGC 6441 in 
contradiction to all scenarios for explaining the upward sloping 
HBs in these globular clusters. Most likely, these large gravities 
can be attributed to problems with the background subtraction 
from the spectra. 

Our results for the hot blue tail stars remain problemat- 
ical. The low surface gravities and masses found for these 
stars might be due to problems with the background subtrac- 
tion or possibly to inadequacies in the stellar atmospheres used 
in the analysis. Further work is needed to clarify this point. 
Unfortunately the atmospheric parameters for the hot blue tail 
stars do not permit a stringent test of the helium pollution sce- 
nario. In particular, the present data cannot test the prediction 
for a greater gravity offset from the canonical ZAHB in the hot 
blue tail stars compared to the cool blue tail stars, as one would 
expect if the helium abundance increased along the blue tail. 

The present study illustrates the difficulties encountered 
when doing spectroscopy in crowded fields. From our experi- 
ence we conclude that the problem of the hot stars in NGC 6388 



cannot be solved by means of current ground-based observa- 
tions. One would need either optical spectra at a much better 
spatial resolution (around (X'l) or multi-colour photometry in- 
cluding the near and far UV. Spectroscopy at such high spatial 
resolution is currently available only at near-lR wavelengths, 
which would worsen the contrast between the hot and cool 
stars and, in addition, would allow only the use of hydrogen 
lines, for which the model atmospheres have not been tested^. 
Multi-colour photometry might provide an estimate of the ef- 
fective temperatures and luminosities of these stars, but as they 
are hot, the satellite UV is definitely required. 
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